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CN ■ Abstract 



Volume conserving surface models (VCS) without deposition and evaporation, as flux free ideal molecular- 
beam epitaxy models, are prototypes to study of the symmetries of conserved dynamics. In this work we study 
two similar VCS models with conserved noise, which differ from each other by the axial symmetry of their 
dynamic hopping rules. We use a coarse-grained approach to analyze the models and show how to determine 
the coefficients of their related continuous stochastic differential equation (SDE). The employed method 
makes use of small translations in a test space, where the test space, which considers the implicit symmetries 
^ \ of the system, contains the stationary probability density function (spdf). In case of the symmetric model 
we calculate all the coarse-grained coefficients of the related conserved Kardar-Parisi-Zhang (KPZ) equation. 
Further we determine the implicit relation between the coefficients and the intensity of conserved noise. With 
respect to the symmetric model, the asymmetric model adds new terms which have to be analyzed, first of 
all the diffusion term, whose coarse-grained coefficient can be determined by the same method. In contrast 
to other methods, the used formalism allows to calculate all coefficients of the SDE theoretically and within 
limits numerically. Above all, the used approach connects the coefficients of the SDE with the spdf and 
O ■ hence gives them a precise physical meaning. 
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1 Introduction 

. The roughening properties of nonequilibrium surface systems with well-defined power-law behavior can be 
classified into universality classes by a unique set of exponents, which determine their dynamical scaling [T]. 
*n , Discrete models and continuous equations within the same universality class share the same scaling exponents, 
00 ' which can be shown by numerical and analytical methods. In order to determine these exponents, many 
discrete models and continuous stochastic differential equations (SDE) have been studied using simulations, 
symmetry analyses, dynamical renormalization group theory, or numerical integration [2]. 

Apart from the calculation of the scaling exponents, there is a particular interest to determine the 
coefficients of the continuous stochastic differential equations, because models and equations within the same 
class share not only the same scaling exponents, but also the same linearities and (or) nonlincarities [3]. In 
this context, models belonging the Kardar-Parisi-Zhang (KPZ) or to the conserved KPZ universality class 
are of special importance, since for these models it is possible to determine the nonlincarities via Monte Carlo 
simulations using the interface tilting method. Another method, which is not limited to the calculation of 
the KPZ-nonlinearities, was introduced by Vvedensky et al. [3]. He starts his calculation from a discrete 
SDE (or discrete Langevin equation) and determines the continuous counterpart of the discrete SDE using 
an analytical approach. The approach is based on the regularization of the Heaviside 0-function and on the 
coarse grain approximation using a lattice constant which is tending to zero [5]-[8]. 

Recently, we introduced a new, different coarse-grained approach based on generalized function or distri- 
bution theory [5]. We showed, that using our approach, it is possible to calculate not only nonlinear, but all 
coefficients of the stochastic differential equation for a given discrete model of the KPZ universality class. 
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In this work we use our formalism to show, how to determine the coarse-grained coefficients of two volume 
conserving surface models. 

There are two well-defined groups of models and equations distinguish by their noise, which is either non- 
conservative or conservative. Growth models and equations which describe deposition or evaporation process 
are included in the first group. In contrast, volume conserving surface (VCS) models without deposition or 
evaporation are included in second group. The models studied in this work are VCS models and consequently 
have conservative noise. Volume conserving processes are defined as physical processes which occur on the 
surface of a solid and preserve the total volume enclosed by the surface. These processes describe the 
movement of a particle from a site of the surface to another, and exclude particle deposition or evaporation. 
The conserved noise is assumed to be Gaussian distributed and uncorrelated. It has an expectation value of 
zero and the correlation is 

(r)(x, t) t}(&, t')) = -2 Q V 2 5{x - x") 5{t - t') , (1) 

with x € ~BL d and where the conserved noise intensity Q is proportional to the temperature of the system. 
The conserved KPZ (cKPZ) equation [h = h(x,t), v < 0, and A < 0] 



dh = 
dt 



+ V(x,t), (2) 



proposed by Sun et al. [TU] describes the first continuous VCS process with conserved noise. The expression 
in the brackets of eq.@ we call KPZ kernel of the equation. In the first molecular-beam epitaxy (MBE) 
models, eq. © reappears with an additive constant F that takes the deposition flux into account. This 
extended equation is called Villain-Lai-Das Sarma (VLD) equation, named by the authors of the article it 
was first time mentioned [HUH]. In contrast to the noise term of the cKPZ equation, the noise term of the 
VLD equation describes non-conservative noise [51113] . It has an expectation value of zero and the correlation 
(rj(x, t) r](x' , t')) = —2 D 8(x — x') 6(t — t'), where the noise intensity D is proportional to the flux F. Although 
the cKPZ and VLD equation have the same KPZ kernel, the equations show distinct macroscopic properties, 
which have been identified and studied in various scientific works [14] . Since both equation have the same 
kernel, these distinct properties can be attributed the different nature of their noise. 

A VCS model can be described by the continuity differential equation or conservation law 

9h „ . 

W + V. J = ,, (3) 

where j = t) is the surface diffusion current and r\ = r](x,t) is the conserved noise. The form of the 
current j is determined by the symmetries of the system. The current cannot depend explicitly on h, since 
this would break the invariance under constant height translations. It is expected that under nonequilibrium 
conditions the total current j has two terms, i.e. j = — V/i nc + j ne - Here fi no is the nonequilibrium chemical 
potential and j no is the nonequilibrium current. From symmetry consideration it is expected that fi nc is a 
function of V 2 h and V/i. 

The first model studied in this work has both of these terms. Its dynamic rules are axial symmetric, 
which means, that a randomly chosen particle hops independently from the actual surface configuration with 
the same probability to left as it hops to the right. Its KPZ kernel [see the eq. ©] is proportional to fi ne . 
Furthermore it is expected that j ne is an odd function of V/i. 

The second model studied in this work has both types of current contribution, V/i no and j no . Its dynamic 
rules are not axial symmetric, which means, that it depends on the surface configuration whether a randomly 
chosen particle hops to the left or to the right. In order to study their models, the majority of theoretical 
studies employ expansions of j in powers of V/i, its derivatives, and their combinations taking into account the 
symmetries of the system. The most relevant term of these expansions is determined by the long wavelength 
limit after renormalization |15j . 

The main goal of this work is to introduce an approach that allows to determine the coarse-grained 
coefficients of the continuous SDE. Additionally, we show how to use the approach to calculate coarse- 
grained coefficients related to two different discrete VCS models, one with symmetric hopping rules and the 
other with asymmetric hopping rules. The paper is organized as follows. In Section [2] we developed some 
basic concepts and the theoretical formalism which is used to study the chosen models. In this context we 
give an overview on the conserved dynamics and conserved noise. In Section [3] we introduce a theoretical 
approach that allows to determine the coarse-grained coefficients using generalized functions or distributions 
theory and proper test functions. We show how to obtain these distributions and the test functions. We 
discuss how small changes of the test functions together with a Taylor expansion of the distribution can be 
used to calculate the coarse-grained observables by apply the expanded distributions to test functions. In 
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Sections@]and[5]we apply the formalism on the chosen discrete VCS models. We obtain analytical expressions 
for the coefficients in terms of test functions. In Appendix A we show how to calculate the coarse-grained 
coefficients using the standard coarse-grained approach via the regularization of 0-functions. In Appendix B 
we evaluate the results of Section H and [5] numerically using Monte Carlo simulations. Section |B] finishes the 
paper with our conclusions. 



2 Volume conserving models 

By definition, volume conserving surface models keep the total number of particles enclosed by the surface 
constant. The movement of the particles can only increase (or decrease) the height of a surface site at the 
expense of decreasing (or increasing) the height of another surface site. Thus, the average velocity of the 
surface growth is zero and expectation value of the surface height is constant. We assume a discrete process 
that takes place on a square lattice with unit cell of side a. Further it is assumed, that a randomly chosen 
particle can move only one unit in a unit time. The surface configuration H of the system is determined by 
set of heights {hj} which correspond to the columns j = 1, . . . , N. The transition rate W(H, H ) of a VCS 
model, which describes the change between two consecutive surface configurations H and H' in the lapse r, 
is 

i N 

W{H,H') = -J2 H A(h' k ,h k - a) A{h' k+1 ,h k+1 + a) 



T 

fe=l 



+ u k+1 A(h' k+1 ,h k+1 -a)A(h' k ,h k + a)] 

x J] A(h'phj) , (4) 

where u> k and 0J k+1 are the hopping rates from column k to column fc + 1, and from column k + 1 to column 
k, respectively. Here A(x, y) is equal to 1 if x = y and equal to in all other cases. The hopping rates depend 
on the height differences between the chosen column and its nearest-neighbor (NN) columns. Explicitly, ur^ 
is function of /ij+i — hj and hj-i — hj. The first transition moment is 

K ? = - M W{H, H') = - r - W 7 - w + + , (5) 

H' 

and the second transition moment, written in an equivalent form to the expression used by ref. |161 eq. 6], is 

K S ] = ~ - hj) W(H, H>) 

ii' 
a 2 r 

= [{ut + + uj+_ x + w i+1 ) (5 i+ ij - 2 S itj + Si-ij) 

+ K + - ^7 - "i-i + wr +i) - s i-hi) > ( 6 ) 

where 6i t j is the Kronecker delta. Note that the first term in brackets contains the symmetric contributions 
and the second term in brackets contains the asymmetric contributions. If the hopping rates away from a 
chosen column i are equal, i.e. = w~, and the hopping rates towards the column i are also equal i.e. 
w^lj = then eq. (0]) has only symmetric contributions. This means, that in this case the jumping rates 

and the second transition moment are invariant to the exchange of hi— n <^- hi +n (with n = 1,2) around to 
column i. 

Starting from the Master equation of the transition probability given by eq. ^ and performing a Kramers- 
Moyal expansion [TTHT^] . it is possible to derive an associated discrete Langevin equation 

^t=K? +„(,), (7) 
with j = 1, . . . , N . The expectation value of the noise rjj is zero and the correlation is 

( Vi (t)r 1j (t>))=Kg ) 5(t-t>). (8) 

Jung and Kim showed |16| that the continuous SDE can be achieved from its discrete counterpart [eq. ([7])] 
following the approach introduced by Vvedensky et al. @J. In Appendix A we use Vvedensky's approach to 
calculate the coarse-grained coefficients of the models which we discuss in the Sections (jj) and (0 in terms 
of the regularizing coefficients. 
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3 A bridge from discrete to continuous processes of the same uni- 
versality class 

According to previous section, the particle hopping rate oj^ from the column j — > j ± 1 is a function of 

the integer set {cj± = (hj±i — hj)/a). Thus, the transition moments Kj and K^' [eqs. ([5]) and (jDJ, 
respectively] depend on the hopping rates and the hopping rates are functions of the height differences 
between the NN columns and the selected column. The identity Uj~ = uj is a fundamental property of the 
studied system which is indispensable to extract all the coarse-grained properties derived from the transition 
moments. Using the transition rates ur^ we define a generalized function W(g) with g G M. 2 . We assume, 
that this generalized function matches with the hopping rates W(f? = (Tj) = oj^ , if g = a-j with (Tj G Z 2 . 

In order to calculate the statistical observables of a discrete processes it is necessary to define test functions 
tp on which the generalized functions of g can be applied. These test functions ip have to have different 
properties for restricted and unrestricted processes. Restricted processes require that ip G T>(1Z 2 ), where D is 
the test space of C°°-functions with compact support and 1Z 2 = e— neighborhood^ 2 ). Unrestricted processes 
require that ip G <S(R 2 ), where S is the test space of C°°-functions that decay and have derivatives of all orders 
that vanish faster than any power of g~ l (a = 1, 2). In the present work we define the test functions on the 
base of the surface configuration. For each column j we have two random integer variables that define a surface 
configuration vector (Tj = (crj_, <7j + ) G Z 2 . The chance to find a column j with the surface configuration 
(Tj at the time t is given by the time dependent probability density function P((Tj,t). One can show via 
Monte Carlo simulations, that the time dependent probability density function converges rather fast. Hence, 
for f > 1 one can interpret P((Tj,t) as a time independent, stationary probability density function P st ((Tj) 
(spdf). The set of surface configurations in the steady state is Z 2 = {(Tj G Z 2 / P st (<Tj) ^ 0, Vj = 1, . . . , N} . 
In the present work the steady state configuration space Z 2 is a 2-dimensional lattice within a compact set 
1Z 2 C (R 2 \ Z 2 ) U Z 2 . We define the test function p(g) as a real-valued function, that matches with the 
discrete spdf ip(g = a-j) = P st ((Tj), if g — a-j with a-j G Z 2 . We used here the notation on distributions that 
was introduced by Schwartz [2D] . 

The application of a distribution / G V (dual space of T>) to test function p G T> is defined by 

(f,<p)= [ /(f?Me)dv e . (9) 

Js. 2 

Note that (/, <p>) is the "expectation" value of / using the test function tp as real-valued analytic "representa- 
tion" of the spdf P st . Thus, the test function is normed, i.e. (1, tp) = 1 . The translation T a of a distribution 
/, denoted T a f or simply f a , extends the definition given by eq. ([DJ 

{T a f,<p) = (f,T- a <p), (10) 

where the translation operator is defined by T x : y i-> y — x if y,x G M 2 [21 . As mentioned above, 
we assume that the test function tp takes fixed values in the discrete lattice Z 2 given by the spdf P s t, i.e. 
{5 a , p>) = p(cr) = -Pst(c) for all cr G Z 2 , where <5 CT = Ta-8 and S is the Dirac distribution. 

Applying a translation T u to a point g G supp(^j), the test function transforms as <p — > T u tp if (g — u) G 
supp(p). Furthermore w = (W, tp), the expectation value of W, changes as w(0) — > w(u) with 

w(u) = {W,T ui p)= [ W(g) <p{q - u) dv e . (11) 
Js. 2 

For small translations, the Taylor expansion of tp(g — u) around u = is 

p(g -u) = p(g) - u a d a tp\ u=Q + \ u a u p d 2 ap tp\ u=Q + 0(3) . (12) 

Here repeated subscripts imply sums, u a is the a-th component of it, d a =d/dg ai and (J 2 ^ =d 2 /(dg a dgp), 
with a ,/3 = 1,2. Since the test function tp is known only at points of the lattice -Z 2 , its derivatives can not 
be calculated explicitly. In contrast, the distribution W is derivable in all points. Since the test function has 
either compact support or decreases rapidly one can take advantage of the following identity 

(W,^..^) = (-l) n (^W^). (13) 
Using eqs. (fTTHI!?]) . the observable is transformed according to 

w(u) = (w, p) + (d a W, ip) u a + %(a£ p W, p) u a u p + 0(3) . (14) 
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In a previous work [5] we used the mentioned translations to determine the coarse-grained coefficients 
for a restricted and an unrestricted discrete KPZ-type model: the restricted solid-on-solid model and the 
ballistic deposition model, respectively. We derive their coefficients from the transformed average velocity 
of the interface v(u) = (K^- 1 ' ,T u ip) where the first transition moment is the drift of corresponding 

Langevin equation. The VCS models studied in the following sections have no restrictions, although the 
theory allows the treatment of models with restrictions, as shown in previous studies. In contrast to this 
former work, the presently studied models allow to calculate analytically another observable quantity which 
is the nonconserved noise intensity. This is possible by intervening in the equivalent continuous process, in 
this case via the application of small changes to D = {K^ 2 \ <p), where is the second transition moment 
corresponding to discrete process and D is the noise intensity corresponding to equivalent continuous process. 



4 Symmetric hopping model (Krug Model) 

At first we study a model with symmetric hopping rate suggested by Krug [14] and studied by Jung and 
Kim [16] . The hopping rules are as follows: 

Chose column j randomly, if the surface configuration satisfies the condition hj±i > hj, the surface 
remains unchanged and a new column is chosen. If the condition hj±i > hj is not fulfilled, the particle on 
the surface of the column j moves with an equal probability either to the left or to the right neighboring 
column. From the rules one can see, that a particle is mobile even if the surface is totally flat. This is necessary 
in order to generate a non-equilibrium contribution since there is no particle deposition. Introducing the 
hopping rates w^, from the chosen column k to right or left NN column, the first moment of the transition 
probability between two surface configurations is given by cq. (JSJ where the hopping rates from column j are 



The first moment of transition rate is 



±Q(h j -h j - 1 )Q(h j -h j+1 ). (15) 



K^=a- 2 ^Z s {g_, Q+ ) (16) 

J Q ± ={hj±i—hj)/a 

where A^ 2 ) is the second variation [i.e. A^Fj = Fj + i — 2 Fj + Fj_i], with the property that a _2 A^ 2 ' — > 
d 2 /dx 2 when a — > 0. Here, the kernel Z s : M 2 — > M. is the generalized function defined by 

Z s ( Q _, e+ ) = ^e(-g_)Q(- Q+ ). (17) 

with a = a 3 /r. Notice that Z s is antisymmetric, i.e. Z s (x, y) = Z s (y,x). The real-valued generalized 
function Z s can be achieved from its discrete counterpart by regularization techniques. Via the coarse 
approximation a stochastic differential equation for h = h(x, i) 

^=V 2 Z s +r,(x,t) , (18) 

can be obtained, where the KPZ kernel is 

Z s = vV 2 h+ ^ \\7h\ 2 . (19) 

Here and below, we use the following notation for the partial derivative V = d/dx. Equation IT51 is a 
conserved Kardar-Parisi-Zhang equation with conservative noise [10] whose expectation value is zero and 
with the correlation function given by eq. ([T]) . In the Appendix A we find this equation using coarse-grained 
techniques via regularization by an alternative way to reference |16j . 

The coefficients v and A of KPZ kernel in terms of the regularizing constants arc given in eqs. (|39[) and 
respectively. 

Instead of using coarse-grained techniques via regularization, the coefficients v and A can be obtained 
applying the distribution Z s [eq. (fTT)) ] on translated test functions T u (p. 

We introduce the transformations g ± — > g ± + u ± and we assume that gradient and Laplacian transform 
with u_j — u +1 = 2 s and it_ 1 + u +1 = 2 r , respectively. Considering a translation u = (r + s,r — s) and 
taking into account that tp is symmetric under exchange g_ -f-> g + one calculates 

2 „3\ 



z s (r, s) = (Z S , T u <p) = uj a + v r + \ A s 2 + 0(r 2 , s 3 ) , (20) 
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Figure 1: The figure shows the test space. In the test space one sees the test function of the symmetric 
hopping model ip(g,,Q_) (semi-transparent) and the discrete probability density function P st (cr+,<7_) 
(black bullets). One sees, that the test function matches with the discrete probability density function 
for q = a. 



where 



W = (Z a ,(p) , 

v = 2 (d e Z B , ip) , 

A = 2((d 2 u -dl)Z s>V ) , 



(21) 



with £ = k — ±1 and I ^ k, where repeated subscripts do not imply sums. Here ip(x, y) is the (symmetric) 
test function. For discrete values x, y G Z the test function ip(x, y) takes the same value as the probability 
density function P st (x,y). Taking in to account that 



we obtain 



CJ 



d £ Z s = 
d'ipZa = 



—\a8{Qi) <d{-Qk) , 
^aS(ge) S(g k ) , 
-\a8'{Qt)@(-Q k ) , 



<p(x,y) dxdy , 



= -a <f(0,y) dy 



¥.(0,0) 



(0,3/) 



dy 



(22) 
(23) 
(24) 



Straightforward one can show, that the noise intensity of correlation function [eq. ([!])] is Q = a 2 z s (u) = 
a 2 wq + 0(1) for all translation u with |u| <C 1. When a — > the noise intensity is 



a 2 wq 



(25) 



Please note, that using our approach we have obtained the three coarse-grained coefficients which characterize 
the conserved KPZ equation. The algebraic signs of these results are in agreement with the algebraic signs 
obtained using the coarse-grained approach via regularization |16| . 
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5 Asymmetric hopping model 



In this section we study a model with asymmetric hopping rate introduced by Jung and Kim |16j . In 
contrast to the previous studied model where a chosen particle moves, independently from the local surface 
configuration with the same probability to the left or to the right side, here a chosen particle moves to the 
NN column which has the lower height. If the heights of the NN columns are both lower than the height of 
the selected column, the particle hops randomly to one of them. Thus the symmetry of the Krug model is 
broken. The hopping rates for this model are 



<4 



= § Q(h k - hh-i) Q(h k - h k+1 ) [1 - sgn{h k±1 - h kTl )] , (26) 



where sgn(z) = 0(z) — 0(— z) is the sign function. The first and second term are the symmetric and 
antisymmetric contributions to hopping rate under the exchange hj-i <H> hj + \ respectively. In contrast 
to the hopping rate w^ 1 , the hopping rate is asymmetric under the exchange hj_ n <-» hj +n (n = 1,2). 
The first moment of the transition rate is 



^A^Z s (g_,g + ) + j-A^Z a (g_,g + ) 



(27) 



where Z s is the symmetric contribution known from the Krug model given by eq. (|17p and is the first 
variation between the NN columns [i.e. A^Fj = Fj + \ — Fj-i] with the property (2a) -1 A2 — > d/dx when 
a — > 0. Just like in the case of the Krug model, the kernel Z a : M 2 — > M is interpreted as a generalized 
function defined by 

= -PQ(-Q+)Q(-Q-) sga(g + -g_) . (28) 

with /3 = a 2 jr. Notice that Z a is antisymmetric, i.e. Z a (x,y) = —Z a (y,x). 

By regularization techniques and the coarse approximation a stochastic differential equation can be 
obtained. In this case h = h(x,t) is 

dh 

— = V 2 Z s + VZ a + V {x,t), (29) 



v Vh + 1 V 6 h, (30) 



where the antisymmetric kernel is 

Z 

and Z s is the KPZ kernel cq. (IT91 with coefficients given by cqs. (|2"2"|) -(| 

We perform a translation u = (s, — s) on the surface configuration space with s«l ; taking into account 
that tp is symmetric in its variables 

z a (s) = (Z a , T uV >) = i/o s + 7 s 3 + 0(5) , (31) 

where z a (s) = and 

= ((dx - d y )Z a , tp) , 
7 = ({d x -d v fZ a ^) . (32) 

Eq. (|31[) contains only odd powers in s since Z & is antisymmetric and contributes only to their odd order 
derivatives which are symmetrical. Below we show only the first term, although the calculations can be 
extended to higher order. Taking into account that 

(d x - d v )Z & (x, y) = [6(x)Q(-y) - S(y)Q(-x)] sgn(x - y) 
-A@(-x)Q(-y)5(x-y), 



wc obtain 



f° 

v = 2p [<p(0, x) - 2 (p(x, x)} dx . (33) 
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Figure 2: The figure shows the test space. In the test space one sees the test function of the asymmetric 
hopping model ip(g,,Q_) (semi-transparent) and the discrete probability density function P st (cr+,<7_) 
(black bullets). One sees, that the test function matches with the discrete probability density function 
for q = a. 

6 Conclusions 

In this work we describe a new approach to calculate coefficients of continuous SDE, which is based on well- 
established concepts of the theory of generalized functions. We apply this approach and show how to calculate 
the coarse-grained coefficients of continuous SDE from two discrete VCS models. For reasons of clarity and 
comprehcnsibility both models are (l-l-l)-dimensional particle models, which differ from each other by the 
axial symmetry of their dynamic hopping rules. For the calculation we take advantage of symmetries in the 
transition moments of the models, which describe the particle hopping rates from one randomly selected 
column to its next neighbor columns. The hopping rates depend only on the height difference between 
the selected and the next neighbor columns. We extend these discrete hopping rates (dynamic evolution 
rules) to the continuum by interpreting them as generalized functions or distributions. Finally we apply the 
generalized functions to test functions, where the test functions are calculated on the base of the spdf of the 
surface configurations. 

The introduced approach has the advantage to unveil the symmetry properties implied by a particular 
model and to determine the coarse-grained observables directly from the evolution rules. Although it is 
impossible to use the approach to determine directly the SDE, it can be used to calculate the coefficients of 
this equation employing a test function which is based on the spdf. In this way the approach gives a very 
precise physical meaning to these coefficients, by relating them to the spdf. 

In contrast using the classical coarse grained approach the obtained coefficients of the SDE are dependent 
of the regularizing coefficients, thus are dependent on the regularizing function and can not be univocally, 
numerically determined. Using our approach we give a precise meaning to these regularizing coefficients and 
determine them univocally relating them to the coarse-grained coefficients of both formalisms. 

We believe that this formalism can be extended to conserved dynamic models with nonconscrved noise, 
such as Wolf- Villain [TT] and Das Sarma-Tamborenea [12] models of MBE. But, the crossover between 
universality classes, the change in the continuous SDE and therefore in their coefficients, suggest that the 
notion of configuration space introduced here must be revised for these models. 

Appendix A 

In this appendix we show how to obtain the continuous SDE starting from the discrete volume conserv- 
ing models and using the coarse-grained approximation via regularization of Heaviside-functions. We use 
a slightly different procedure than other authors, but obtain the same coefficients and equations. The ex- 
pressions obtained here are related with the expressions obtained in Sections HI and [5] by the coarse-grained 
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approximation via generalized functions introduced in Section [3J 

Using the regularization procedure, the Heaviside function Q(x) can be replaced by a smooth real- valued 
function e (x), which satisfies 9 £ (n) — > Q(n) when e — > + for all neZ. There have been several proposals 
to represent for a shifted analytic function, including the following 



erf 



erff- 



ds 



with e > 0, introduced in refs. [221125] that result well. The kernels Z s and Z a [eqs. ([XT]) and ([2"8")l . respectively] 
can be regularized using the e-theta function 8 £l i.e. 



Z s {g_,g + ) = la0 s (-Q-)0e{-Q+) 

z a { Q _,e + ) = -pe e {-Q_)e e {- e+ )[e e { e+ - e _) -6 e { e _ - q,)] 



with a = a 3 /t and j5 = a 2 jr. Expanding the e-theta in Taylor series around x = 



fe=0 



gives the following expansions (superscripts are omitted hereafter) 



Z a (g_,g + ) 



+ 



Al - A A 1 (g_ + g + ) 

Q 2 _ - ^1Q- Q+ + Q 



(2A A 2 - Aj) 
2(7 + 1) 



2 ' 

1 



0(3) 



where 



Evaluating eqs. ([35 



Z s (g_,g + )\ 
Z a (g_,g + ) 



2A 2 A 1 ( Q+ -g_)+O(3) 
Aj 

7 2A A 2 ■ 

(i) 



= W0 + ^r + ^^ 7) + O(3) 



e ± =('ij±i-''i)/o 



^o^ (2) +0(3) , 



where the constants are 



(34) 



(35) 



(36) 



(37) 



UJQ = 

V = 

A = 

= 



2r 



-.4 



! 



£ T 



n 



(2A A 2 - Af) 



2r 

AqAi , 

T 



and the linear and quadratic terms are 

(i) _ hj+i - 2hj + 



(2) 



h 



iV 



(7) 



2a 

Pj+i - fej) 2 - 2 7 - hjjjhj-i - fej) + (fej-i - 
2a 2 ( 7 + l) 



(38) 
(39) 
(40) 
(41) 



(42) 



with < 7 < 1 the discretization parameter of the discretized nonlinear term [7]. The usual choice 7 = 1 is 
called standard or post-point discretization. It depends only on the height of the NN columns and thus the 
error of approximating (V/1) 2 is minimized. In contrast, the choice 7 = 0, called antistandard or prepoint 
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discretization. It corresponds to the arithmetic mean of the squared slopes around the interface sites. In 
this work, 7 depends on the coefficients of regularization through eq. ([55]). as Ao > then A2 < 0, that 
coincides with the result obtained in ref. |TS] . In addition we found if A\ > 0, then v < 0, A < 0, and v$ < 0. 
Expanding eqs. (|42[) around x = ja, the discretized terms and their limits when a — >• are 

= V 2 h + — V 4 h a 2 + 0(4) — ► V 2 h , 

J 12 

L {2) = Wh+ ^-V 3 ha 2 + 0(4) — ► Wh , (43) 

J 16 

iV, (7) = (V/i) 2 + -(\^-\w 2 hf a 2 + 0(4) — > (V/.) 2 . 
J 4 \1 + 7/ 

Notice that the limit of iV. does not depend on the discretization parameter 7, as shown in ref. [7J. Take 

into account that — Z s — > Z s and —Z a — > Z a when a — > 0. Applying these limits to eqs. (TTB| and (|27[l we 
obtain eqs. (|18[) and (f2l)]). respectively. A dimensional analysis of these last continuous differential equations 
shows that the their coefficients given by eqs. ([38|) - (j4Tj) have the correct dimensions. In order to obtain the 
continuous limit of cq. © we take a — > 0. The following limits are calculated for eq.©: 

a- s (&t+ij-2&i ti + &i-xj) —> V 2 S{x~x r ), 
a 5 T-\ui + ur +u} +_ 1 + ul - +1 ) 4 Q. 

Using these limits we conclude, that K^j — > K^ 2 \x — x') — —2Q\7 2 S(x — x'), and the continuous noise 
correlation function is given by eq. ((T|). Finally, the noise intensity is 

Q = f^. (44) 

Appendix B 

In this appendix we give the results of Monte Carlo simulations performed for the studied volume conserving 
models. We calculate the coarse-grained coefficients via theory of generalized functions. Both studied 
models evolve over time unit r = 1 in a square lattice with unit cell size a = 1, with a total system size N 
and use periodic boundary conditions. On the discrete lattice we calculate the stationary joint probability 
distribution -P s t (c_ , er + ) of the surface configurations with a ± = hj±i — hj in order to evaluate the coarse- 
grained coefficients (see Table [I) . The integrals of the test function ip or their derivatives are evaluated by 
conventional finite difference methods, based on P st data. With the coarse-grained coefficients of Table [1] we 



Model 


^0 


V 


A 


^0 


Q 


SHM 


0.180 


-0.116 


-0.034 




0.180 


AHM 


0.250 


-0.193 


-0.158 


-0.615 


0.250 


Eqs. 


(E2D 


(123]) 


dSD 


022) 


422), 025) 



Table 1: Coarse-grained coeficicnts for the Symmetric hopping model (SHM) and Asymmetric hopping 
model (AHM), calculated according to equations given in the 4-th row. 

calculate directly the coefficients of regularization and discretization parameter (see Table ^ . 



Model 


A 


Ax 


A 2 


7 


SHM 


0.600 


0.139 


-0.040 


0.406 


AHM 


0.707 


0.272 


-0.171 


0.307 



Table 2: Regularizing coefficients of the e-theta expantion and the discretization parameter 7 calculated 
by eqs. (|5^)) - (j4T|) and eq. (f3"6")l . respectively, using the data of Table [TJ 
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